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Abstract 

We show that classical thermodynamics has a formulation in terms of Hamilton- 
Jacobi theory, analogous to mechanics. Even though the thermodynamic vari- 
ables come in conjugate pairs such as pressure/ volume or temperature/entropy, 
the phase space is odd-dimensional. For a system with n thermodynamic de- 
grees of freedom it is 2n -I- 1-dimensional. The equations of state of a substance 
pick out an n-dimensional submanifold. A family of substances whose equa- 
tions of state depend on n parameters define a hypersurface of co-dimension 
one. This can be described by the vanishing of a function which plays the role 
of a Hamiltonian. The ordinary differential equations (characteristic equations) 
defined by this function describe a dynamical system on the hypersurface. Its 
orbits can be used to reconstruct the equations of state. The 'time' variable 
associated to this dynamics is related to, but is not identical to, entropy. After 
developing this formalism on well-grounded systems such as the van der Waals 
gases and the Curie- Weiss magnets, we derive a Hamilton-Jacobi equation for 
black hole thermodynamics in General Relativity. The cosmological constant 
appears as a constant of integration in this picture. 
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1 Introduction 

General Relativity is surprisingly similar in many ways to Thermodynamics. It 
predicts an irreversibility in the dynamics of black holes: the area of the horizon 
must increase. Since we cannot see the interior of a black hole, they carry an 
entropy which must be proportional to this area. The analogy suggests that GR 
is the average over some as yet unknown microscopic degrees of freedom. 

The situation is analogous to that with the theory of gases at the end of 
the nineteenth century. It was clear that the kinetic theory of gases gave a 
good derivation of the ideal gas law and through the work of Boltzmann, a 
possible derivation of thermodynamic laws from the first principles of mechanics. 
Yet, it was not known then what the value of Boltzmann's constant was, or 
equivalently, how big molecules are. Yet , van der Waals was able to deduce a 
simple equation for the departure of a gas from ideal by simple considerations 
of molecular dyamics. In a similar way, it might be possible to test microscopic 
theories of gravity by their implications on the thermodynamics of black holes 
or of the universe. We find already an intriguing connection between black hole 
thermodynamics and the cosmological constant. 

This motivation led us to questions about thermodynamics that turned out 
to be of independent interest . In particular, we found a dynamical system on 
the thermodynamic phase space that is analogous to classical mechanics but is 
still quite different. The 'time' variable of this evolution of the thermodynamic 
variables is related to, but is not always the same as, entropy. Although many 
of the ingredients of our theory were known in the nineteenth century, this 
particular form of thermodynamics appears to be still new. We have tried to 
make the paper self-contained by reviewing some of the background material 
on mechanics (section 1), thermodynamics of gases (section 2) and magnets 
(section 3) and contact geometry (section 4). The erudite reader could skip 
ahead to sections 2.3, 3.3, 5 and 6 for the new material. 

1.1 The Three Classical Disciplines 

Three important disciplines of classical physics are Mechanics, Ray Optics and 
Thermodynamics. The work of Hamilton and Jacobi showed that Mechanics 
and Ray optics had a very similar structure: they can both be formulated in 
terms of the Hamilton-Jacobi (or eikonal) equation. We now understand why 
that is so: they are both approximations to wave theories in the limit of short 
wave length. 

Let us invert history and start with the wave equations of optics and quantum 
mechanics [U [5] 



c2 dt^ 




;?2 



(1) 



If we put (f) — e'>^ ,tp — e'^ri these become in the limit A — > 0, ft — s- 0: 
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These are the Eikonal and Hamilton- Jacobi equations. The phase of the wave is 
constant on the wave front ( the surface S = constant) and varies most rapidly 
along the vector ViS* which is normal to the wavefront. The wave can be thought 
of as propagating along the ray that solves the ODE ^ —\7 S. 

1.2 Hamilton's Equations 

Thus VS* is proportional to the wavenumber (momentum). We can identify also 
^ with frequency (energy) up to a constant. A Hamilton- Jacobi equation of 
the form 

H{x,\/S)^— (3) 

is equivalent to the system of ordinary differential equations that determine the 
rays|3j 

dx dH dp dH 

H " 'dp' H " ~~dx' ^ ' 

These in turn are equivalent to Newton's laws of mechanics ( or Snell's law of 
refraction) . 

Position and momentum appear symmetrically in the above equations. In 
fact there is a much bigger symmetry. Any observable F{x,p) generates a 
one-parameter family of co-ordinate transformations through Hamilton's equa- 
tions. These are the canonical transformations generated by that observable. 
Time evolution is simply the canonical transformation generated by the energy 
( Hamiltonian) function. Momentum generates translations, angular momen- 
tum generates rotations etc. Observables come in conjugate pairs generating 
translations on each other [2]. 

1.3 Classical Thermodynamics 

In this paper we will show that thermodynamics is a third facet of this classical 
framework. Gibbs [4j realized that the basic law of thermodynamics, combining 
the first and second laws, has a geometric meaning [3|. The theory of differential 
forms was developed by Pfaff, Caratheodory and others to give his insight an 
elegant mathematical formulation. Chandrashekhar [5] gives a rare treatment 
of this subject in a textbook. An analogy of thermodynamics with mechanics 
and optics is mentioned occasionally O [3 [9l [10] . The most detailed geomet- 
rical formulation is in the review article Ref. [T^]. It is shown there that a 
hypersurface in a contact manifold defines a dynamical system. 

But the physical meaning of this dynamics appears not to be clarified any- 
wher^. What is the meaning of such a hypersurface in thermodynamics? It is 
certainly not the submanifold defined by the equation of state of a material, as 

^This analogy with mechanics was noticed in Ref. [9], but it is not as perfect as implied 
there. In particular, because the thermodynamic phase space is odd dimensional (see be- 
low), the correct formalism uses contact geometry and Legendre transformation rather than 
symplectic geometry and Poisson brackets. 
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it has the wrong dimension. (More precisely, it is not a Lagrangian submani- 
fold. See below). What is the Hamilton- Jacobi equation of thermodynamics for 
some familiar systems like gases or magnets? What is the meaning of the 'time' 
variable in this dynamics? What arc the implications for the thermodynamics 
of black holes? These are the questions we will begin to address in this paper. 
After understanding these down to earth system we will turn to the dynamics 
of the simplest black hole. 

2 Thermodynamics of Gases 

To be concrete let us consider the most familiar thermodynamic system, a 
gas with a fixed number of molecules. It has five thermodynamic variables 
U, T, S, P, V: the internal energy, temperature, entropy, pressure and volume. 
S,U,V are extensive variables] i.e., arc proportional to the size of the system. 
P, T are intensive: they do not scale with the size of the system. Two subsystems 
are in equilibrium if they have equal values for the intensive variables. 

Extensive and intensive variables come in conjugate pairs such as pres- 
sure/volume or temperature/entropy. Yet there are always an odd number 
of variables. This can be seen in the first law of thermodynamics which is a 
condition on the infinitesimal variations of these quantities: 

a = dU- TdS + PdV = 0. (5) 

The conjugate pairs appear together as PdV or TdS,except for internal energy 
which has no conjugate. 

An equivalent form would be 

dS - T-^dU - PdV = (6) 

where P = ^ . Indeed any variable can be chosen to be in the privileged position 
as the 'fundamental variable' in the first term, with unity as the conjugate. 
This symmetry resolves the apparent contradiction of having an odd variables, 
while at the same that variables come in pairs, a is only determined modulo 
multiplication by a non-zero function: only its kernel (zero set) has a physical 
significance. 

Equivalent formulations of the first law are 

d[U + PV] - VdP - TdS = 0, d[U - TS] + PdV + SdT = (7) 

etc., in which the positions of some of the extensive and intensive variables are 
switched. These transformations are analogous to the canonical transformations 
of classical mechanics and are called Legendre transformations. More on this 
later. 

2.1 The Ideal Gas 

The particular two dimensional surface corresponding to each material is given 
by three equations of state to be satisfied among the variables (T, P, s, u, v). For 
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example, a monatomic ideal gas satisfies the laws 

Pv = T, 



u = —1 
2 



(8) 



where V = nRv is the volume of the gas, U — nRu the internal energy, n is 
the number of moles of the gas in the sample and R w 8.3JK^^mol^^ is the gas 
constant. The remaining relation needed to fix the two dimensional surface of 
the ideal gas is obtained by integrating the differential equation a = 0. This 
gives the entropy E|: 



P 1 

ds — —dv — —du 
T T 







log 



(9) 



This holds up to an additive constant of integration; i.e., independent of 

Conversely, once this 'fundamental relation' between extensive variables is 
given the intensive variables are determined by 



(10) 



1 




P 


(-) 






T ^ 


\dvj 



which are the other equations of state. 
2.2 van der Waals Gases 



A simple model of a non-ideal gas is due to van der Waals [14]: 

T 

> — b 



P 



a 



(11) 



where a, b are constant parameters. We can regard 6 as the volume excluded 
because of the finite size of the molecules. The constant a measures the strength 
of the short range attraction among the molecules. 
Let us now find the entropy: 



ds — 



P 



-dv 



dv 



du 
Y 

du %dv 



(12) 

v-b T ^ ' 

Thus T is an integration factor which must make the last term an exact differ- 
ential. A moments thought gives 

2 r 



(13) 



We fix the constant of integration by noting that the low density limit v —> oo 
must be the ideal gas. Thus, 



log 



{v-b)[u+- 
* v 



(14) 



= nRs, where n is the number of moles of the gas and R the gas constant. We will 
find it useful similarly to define 'specfifor all extensive variables S = nRs, V = nRv, U = 
nRu, M = nRm etc. The gas constant is put in to simplify later formulas. 
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To check: 



P 
f 



(-) 



3a 



21)2 



1 

T 



[ihi 



2 u - 



from which it follows that 



(15) 
(16) 



Now we can eliminate the two parameters a, b from the three equations 
(|11I13I16[) to get an equation relating the thermodynamic variables: 



P 



1 



— V u H — 

T T 2 



rr-3 27 2 - 

1 = —V e 



(17) 



This is the hypersurface in the thermodynamic phase space describing the 
whole family of van der Waals gases. We want to find such an equation for all 
the examples of thermodynamic systems we study. 



2.3 Hamilton-Jacobi equation for van der Waals gases 

From the equation (|17p we can now recover the equations of state of the van 
der Waals gases. In the picture in which s is the thermodynamic potential and 
u, V are the co-ordinates we can put in for their conjugates 



ds 
du 



to get the first order PDE 



ds 



ds 



dv du 



P 
T 



ds_ 

du 



27 



(18) 



(19) 



This is the Hamilton-Jacobi equation of the thermodynamics of van der 
Waals gases. The complete integral of this equation ( in the sense of Ref.[3]) 
will have two undetermined parameters. This is exactly ((1] 



3 Magnetic systems 

For a magnetic system, the analogue of pressure is the applied magnetic field B 
and its conjugate variable is the magnetization M . The law of thermodynamics 
become^ 

a = dU - TdS - BdM = 0. (20) 

The equations of state that determine the two dimensional surface of a mag- 
netic material is known if we know S{U,M) or the Gibbs potential G{B,T). 
We expect that the entropy will decrease with increasing magnetization as the 

•^For simplicity we consider only one component of the magnetic field and magnetization. 
When the magnetic susceptibility is isotropic this is sufficient. 
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system is getting more ordered: its molecules are more aligned. Again by consid- 
ering one mole of the paramagnetic substance, the thermodynamic law becomes 



1 13 
ds — —du — bdm, ^ ^ j^' ('^■'■) 

where, as before, 

S = nRs, M = nRm, U = nRu (22) 
and n is the number of moles of the magnetic material. 

3.1 The Ideal Paramagnet 

The magnetic analogue of an ideal gas is an ideal paramagnet, consisting of 
dipoles that don't interact with each other, only with an applied external field 
[7]. Thus the internal energy of the magnets is independent of m, same as what 
it would have been in the absence of magnetism : 

u = uoiT). (23) 

In these variables b,m,T , temperature 'decouples' from the magnetic variables 
and we can treat the two degrees of freedom independent of each other. 

ds = -6dm + ^diio(r). (24) 

Each material has an equation of state that relates 6 to m . We would expect 
that TO is an odd function of b, which tends to some finite limiting value toq as 
6 — > oo: when all the molecular magnetic moments have been aligned. A simple 
model that captures this idea would be 

to — Too/ {bmo) (25) 

for a constant toq and an odd function /(y). The function f{y) should tend to 
±1 as 2/ — > ±oo. Also, f{y) ~ y + 0{y^) for small y . If 5 is the inverse function 
of /, so that gifiy)) = y, 

6 (26) 

Too ymoj 

The magnetic part of the entropy would then be given by ^ = — 6 so that 



si(to) — — I g{x)dx. (27) 

Thus 







si{m) + so{T) (28) 



where = i^«o(T) 
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For example, the case of spin | magnetic moments 

f{y) = tanhy, g{x) = ^ log (^^^f ) ' " -^x\og ~\ ^og[l-x^]. 

(29) 

But a simpler choice leading to more explicit answers is 

f{y)=y^, gix) = ^=l=, S^{x) = Vn^-l. (30) 

These two functions f{x) and g{x) are similar enough to tanhx and arctanh x 
that the simpler choice already gives us a good idea of what is going on. The 
simplest model for the non-magnetic part of the internal energy would be that 
the molecules carrying the magnetic moment are oscillating around the vertices 
of a lattice. If there are r such oscillatory degrees of freedom, equipartition of 
energy gives 

uo{T) = ^-T. (31) 

Thus 

so(T) = ^logr. (32) 
3.2 Curie- Weiss Theory 

Just as the van der Waals theory gives a simple model for non-ideal gases, 
the Curie- Weiss theory gives a simple model for ferromagnetism. Each spin 
is thought of as interacting with a local magnetic field that is the sum of the 
applied field and a constant multiple of the average of all the other spins. Thus 
each model of ideal paramagnct can be turned into a Curie- Weiss model for 
ferromagnetism. The formula for magnetization is modified to the implicit ( 
self-consistent) formula 



m = mo/ (mo6+ ) 
V T mo J 



(33) 



The parameter Tc ( the Curie temperature) measures the strength of the inter- 
actions of the magnets. When 6 = 0, the solution remains m = as long as 
T < Tc- If T > Tc, there can be spontaneous magnetization even when 6 = 0. 
By writing the equation of state as 

1 / m\ m 



mo \ mo / T mg 
and requiring that 

ds = -hdm + — (35) 

be an exact differential, we get 

1 f Tfl \ ^ / Tfl \ 

u{m,T) = --TA — \+uo{T), s = sA — ]+so{T) (36) 

2 \mo J \mo J 
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where uq and si arc the same functions as for the ideal paramagnct. 

We can now eUminate the parameters mo,Tc to get a relation between 
s, T, b, u, m that characterizes the Curie- Weiss models. First of all, 



1 771/ 

mb = xg{x) - x^T^-, ^ " m^' ^^''^ 



We eliminate Tc using , 



u = uo{T) - iar^Te (38) 



so that 

mb = xg{x)+2^^^^^. (39) 
If we define the functions h and j by 

h{s,{x))=x, j{y) = h{y)g{h{y)) (40) 

we get 

x = h{s- so{T)) (41) 

and 

mb = j{s-so{T)) + 2'^^^. (42) 

Thus the two functions so(T) and g{x) characterize the mean field model, re- 
calling 

uo{T) = [ Tdso{T), si{x) = - [ g{x)dx. (43) 
Jo Jo 

We can get an explicit form of this constraint with the simple models 

Wo(r) = ^T, si(x) = 1 - Vl-a:2 ^ (44) 

g{x) = ^^=, %) = yT^7^^, j{y)=^--y, (45) 

V 1 - ^2 y 

= - - + 2^i^|^ (46) 



And finally 



1 7" "ZX 



3.3 Hamilton- Jacobi Theory for Curie- Weiss Magnets 

It is simplest to use T, m as the thermodynamic variables. Then the law of 
thermodynamics becomes 

dE = -SdT + BdM, E = U-TS (48) 
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where S is the Hehnholtz Free energy. Defining S = nR(^ as before 



-sdT + Tbdm, = ~ I ^ 



T \dm 



The Hamihon-Jacobi equation becomes, with rj = ^ + ^ 

drj 



m drj drj 2r rj 



dT 



= 1. 



(49) 



(50) 



Again we can find the general solution with two arbitrary parameters which will 
give the equations of state. 



4 Contact Geometry 

We can abstract out of the above a mathematical structure that captures the 
essence of thermodynamics. 

Consider a local patch with co-ordinates (9", 9^ • • • , ^"jPi, • • - Pn) of a man- 
ifold with dimension 2n + 1. A contact structure in this patch is given by the 
one-form 

a ^ dq'^ — Pidq^ , ■ ■ ■ i = 1, • • • n (51) 

up to multiplication by a non-zero function. That is, we are to regard a and fa 
(for a non-zero function /) as equivalent. It is the vanishing of the infinitesimal 
variations that defines a contact structure. 
We can think of 

dq° - pdq' (52) 

as the condition for maximizing (or minimizing ) some quantity q^ subject to 
the constraint that some others are held constant; the pi are the Lagrange 
multipliers that enforce the constraints. Such problems arise not only in sta- 
tistical mechanics and thermodynamics, but also in other areas of physics and 
even in economics [TU[ ITT] . 

There are many other co-ordinate systems in which also the contact form 
(up to multiplication by a non-zero scalar) has the above canonical expression. 
From 

- pdq'] = dQ" - PdQ' (53) 



we get the conditions 



" ^ dp, ' dq' ' dq' ~ 



(54) 



^ ' dq^> dq" 

Transformations {q^^^q^yPi) {Q''\Q^,Pi) that satisfy these conditions are 
called Legendre transformations. 
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A contact manifold is a union of co-ordinate patches such that the trans- 
formations among co-ordinates defined at the intersections are Legendre trans- 
formations. In most cases the contact manifold of interest is just i?^""*"^ with 
the standard contact form above. But there are exceptions: a superconducting 
circuit with a Josephson junction has x as the thermodynamic phase 
space[Tn]. Even on R'^, there are contact structures that are not equivalent to 
the standard one; i.e., that need several patches to cover R^. But it is not clear 
yet whether they are of interest physically. In any case, we will be content with 
a local description in this paper. 



4.1 Generating Functions of Legendre Transformations 

An infinitesimal Legendre transformation 

<7°^g° + iF°, q'^q' + tV\ p^^p^ + m, \t\«l (56) 
defines a vector field 

whose components satisfy the infinitesimal version of the above condition: 



dpi 



Pj 



dpi ' dq' 



I 



9g0 



vr 



V, = p. 



(58) 



(59) 



After a little work we can see that all the components are expressible in terms 
of the single function F = pjV^ — y°: 



dF 
dpj 



dF 



dF 
q^ 



Pi 



dF 
dpt 



F. 



(60) 



Thus an infinitesimal Legendre transformations is determined by a single func- 
tion F, called its generating function. 

Conversely, given a function F, define the vector field 



Vf 



p^-F 

Pi ^ 

Opi 



d 

dq" 



dF 



Pi 



dF 
d^ 



_d_ 

dpi 



dF d 
dpi dq' 



The infinitesimal transformation of the contact form 

a — dq" — Pidq^ 
under this vector field can be calculated using 

CyCi = d[ivoi\ + iyda. 



(61) 



(62) 



(63) 
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For us 



-F, da = ~dpi A dq^, ivpda 



dF 



dF 



dq dq^ 



It follows that 



r 



dq' + ^dp,. (64) 



(65) 



Thus the vector field Vp changes the contact form only by an overall multiplica- 
tion by a scalar: it leaves the contact structure unchanged. Such transformations 
are called Legendre transformations. 



4.2 Characteristic Curves of a Function 

Finite Legendre transformations can be constructed by composing such infinites- 
imal transformations; i.e., by determining the integral curves of this vector field. 

Thus, a generating function defines a one-parameter family of Legendre 
transformations, determined by solving the Ordinary Differential Equations: 



dq' 
~dt 



dF 
dpi' 



dpi 
dt 



dF 
dq' 



dF 



df_ 

'dt 



Pi 



dF 
dpi 



F. 



(66) 



These are called the characteristic curves of the Generating Function F. 

These equations are the analogues of Hamilton's equations in mechanics. 
But note that they are not quite the same: there are additional terms in p if 
the generating function depends on q^. Moreover, the value of F is not always 
conserved: 

(67) 



dt ^ 



drf 



But if the initial value of F is zero, it remains zero. 

Conversely, every hypersurface can be thought of as the set of zeros of some 
function F. The characteristic curves of F define a dynamical system on ev- 
ery hypersurface on a contact manifold. These curves are unchanged ( up to 
reparametrization) by a change F — > ^(i^) of the function defining the hyper- 
surface. This dynamical system will play an important role in our approach to 
thermodynamics. 

An important example of a Legendre transformation is one that interchanges 
a co-ordinate with its conjugate: 



Pi = q^,Q^ = -Pi,P2 =P2,f = g^ 



(68) 



The characteristic curves of the function F{p, q) — ^\pf + (q^)^] are circles: the 
above Legendre transformation corresponds to a rotation through | . 
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4.3 Lagrange Brackets 



The commutator of two vector fields satisfying the conditions Cycy. — gya, Cija = 
gjjCK is also an infinitesimal Legendre transformation. If the generating func- 
tions of U and V are F and G, what is the generating function of [?7, F]? A 
straightforward calculation shows that it is given by the Lagrange bracket 

iF r\-r—-—F f^^^^^\ d£dG^d£dG 

^ ' ^ " dqO~ dqO \dp, dq^ ~ dq° Op. J ^ Op., dq' ~ Oq' Op, ^ ' 

These brackets are different from Poisson brackets of mechanics in important 
ways. Although it defines a Lie algebra ( i.e., antisymmetry and Jacobi identity 
are satisfied), 

{F,G} = -{G,F}, {F,{G,H}} + {G,{H,F}} + {H,{F,G}}^Q (70) 
the Leibnitz rule of derivations is not satisfied in general: 

dF 

{F, GH} - G{F, H} - {F, G}H = -GH— ^ 0. (71) 
Thus for example, even a constant can have non-zero brackets: 



r)G 

(l.G)^-^. (72) 



We have analogues of the canonical commutation relations 

{p„q^}^S}, {q\q'} = = {p^,p,} = {q°,q'}, {q°,P^}^-p^■ (73) 

But due to violations of the Leibnitz identity these have to be used carefully. 
For example, 

{q',P^q'} = 0. (74) 



4.4 Lagrangian Submanifold 

A Lagrangian submanifold 1 1 2| of M is a submanifold of maximal dimension 
(i.e., n), all of whose tangent vectors are annihilated by a . As an example, the 
submanifold 

Pi = 0, • • - pn = 0, 9° = constant (75) 

with co-ordinates (9^, • • •</") is Lagrangian. But not every rt-dimensional sub- 
manifold is Lagrangian. For example, the n-dimensional submanifold with co- 
ordinates {pi,q^ ,q'^ ■ ■ ■ q^) for which 

q*^ = constant, q^ = constant, Ps — ■ ■ ■ Pn — (76) 

is not Lagrangian. The independent variables on the submanifold must have 
zero Lagrange bracket: they cannot include a conjugate pair. Thus a Lagrangian 
submanifold is the analogue of a configuration space in mechanics. 
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More generally, a submanifold is determined by a function 

gO = $(g\---,g"). (77) 



The remaining equations determining the submanifold are 

It is straightforward to verify that any infinitesimal variation within this surface 
satisfies dq° — Pidq^ = 0, since 

dl' = ^dq\ (79) 

The trivial case of a Lagrangian submanifold mentioned above corresponds 
to the choice $ = constant. 

The same Lagrangian submanifold can have different descriptions, as we can 
use different canonical co-ordinate systems on the thermodynamics phase space. 
These different descriptions are related by Legendre transformations. 



5 The Geometry of Thermodynamics 

Now we turn to the formulation of thermodynamics in terms of contact geome- 
try. As noted earlier, if a substance has n degrees of freedom, its thermodynamic 
phase space is 2n + 1 dimensional. The first law can be see as defining a contact 
structure on the thermodynamic phase space: 

a = dq° - pidq' = 0. (80) 

We can, as a first step, think of the 'co-ordinates' as extensive variables and 
the 'conjugate momenta' Pi as the intensive variables. However as noted above, 
a Legendre transformation can mix these up. Having a picture that does not 
depend on the choice of co-ordinates in the thermodynamic phase space can be 
very useful in clarifying complicated thermodynamic relations. This is similar to 
the situation with canonical transformations in mechanics which allow momenta 
and co-ordinates to be mixed. 

n the case of a gas, some of these equivalent descriptions, related by Legendre 
transformations, are 

a = dU-TdS + PdV 

= dE + SdT + PdV, E = U-TS 

= dG + SdT - VdP, G = U -TS + PV 

~ dS-^dU-^dV. (81) 

S is the Helmholtz Free Energy which is the convenient quantity to study a gas 
at constant temperature and volume. The Gibhs Free Energy G is useful to 
understand a gas at constant pressure and temperature. 
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5.1 Equations of State=Lagrangian Submanifold 



The equations of state of a substance define a surface of dimension n in the 
thermodynamic phase space. Thus, among the 2n+ 1 thermodynamic variables, 
there must be n + 1 equations of state. But it cannot be any n-dimensional 
submanifold: the equations of state of a substance must define a Lagrangian 
submanifold. This is simply the condition that any infinitesimal change in the 
state of the substance must satisfy the first law. That is, any tangent vector 
field will annihilate a. 

Once the 'fundamental relation' giving = ^{q^ ,■ ■ ■ q") is given, the re- 
maining n equations of state follow by differentiation as above. As an example, 
for a gas, the internal energy as a function u{s, v) of the other extensive variables 
give the other two equations of state: 




Thus each substance has a Lagrangian submanifold corresponding to it. But the 
same Lagrangian can have many descriptions in different canonical co-ordinate 
systems. The advantage of the geometric point of view is that it allows us to 
choose variables according to the convenience of the physical problem. Often, 
by a Legendre transformation we will be able to simplify equations that need 
to be solved. See below for examples. 



5.2 Family of Substances— > Hamiltonian 

A family of substances which have similar equations of state (e.g., van der Waals 
gases with various values of the parameters a, b ) can be described by allowing 
the function $ to also depend on n parameters ai, • • • a„: 

q° = ^q^...q^\ai,---an). (83) 

These parameters could be quantities such as critical pressure and temperature. 

We require the function defining such a family to satisfy a non- degeneracy 
condition: 

det-— ^0. (84) 
oqoa 

The remaining equations determining the submanifold are 

Pi = $i(g\---g"|ai,---a„), = ^ (85) 

Given the non-degeneracy condition, we can eliminate the parameters ai • • • o„ 
from the above n+1 equations to get a single relation among the thermodynamic 
variables: 

F{q^q'---q'',Pi,---Pn) = (86) 

We will see that this function F{q, p) defines a dynamics on this hypersurface; by 
analogy to mechanics, we call this the hamiltonian of this family of substances. 
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There are some important differences from mechanics because the physics is 
determined by the set of zeros of F. A change of the hamiltonian F (t>{P) 
by an invertible function (j> : R ^ R will leave the surface unchanged and the 
physics will be the same as well. However, adding a constant to F will change 
the physics, unlike in mechanics. 

Thus the thermodynamics of a family of materials is described by a contact 
manifold M, [a] and a hyper-surface F{q,p) = of co-dimension one on it. 

Given one member of the family ( e.g., ideal gas or paramagnet) is there a 
mathematical way of guessing the whole family? Unfortunately the answer is 
no. Such a generalization is essentially a postulate about on the macroscopic 
effects of the underlying microscopic system. Some physical input is needed to 
make the correct generalization. The arguments of van der Waals and Curie- 
Weiss were great theoretical leaps at a time when very little was known about 
molecules. Given a microscopic theory, we can in principle derive ^{q\a) as 
an 'effective field theory'. In the absence of a microscopic theory, some mix of 
experimental information and physical intuition is needed. 

Why do we define a family of substances to have exactly n parameters? With 
less than that many parameters, we are not talking of a 'generic' member. The 
non-degeneracy condition above cannot be satisfied and we would get a surface 
of co-dimension greater then one: it would not be defined by the vanishing of 
a single hamiltonian. On the other hand, with more than n parameters, we 
have the opposite problem: they cannot all be independent, as the surface is 
living in an ambient space of dimension only 2n + l. Thus in a sense additional 
parameters are 'irrelevant': they should not be important for the macroscopic 
dynamics of the system. Indeed in can der Waals theoy, exactly two parameters 
a, b give a good description of the non-ideal behavior of many gases. 

5.3 The Hamilton- Jacobi Equation of a Family 

Given the generating function F(q,p) of a family of substances, we can recover 
the equations of state of each member of the family by solving the first order 
partial differential equation 



The complete integral[3j of this equation will depend on n parameters oi, • • • Un- 
The different equations of state of the family members are given materials by 
different choices of ai,---a„. It is a fundamental tenet of Hamilton- Jacobi 
theory that solving a first order PDF is equivalent to solving a system of ODFs 
for the characteristic curves. 




(87) 
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5.4 The Dynamics of Thermodynamics 



Thus, given the generating function F{q,p) of a family of substances, we get a 
dynamics on the hypersurface F{q,p) — 0. We gel|f|[3j 

dq' _ dF dp^ _ dF dF dq° _ OF 

Ht^d^,' ~dt "W'^'W' ~dr^^'df,' ^ ' 

The parameters characterizing the material oi • • • a„ are given by the initial 
conditions of this dynamics. If we eliminate the 'time' variable we can get 
the complete integral of the Hamilton-Jacobi equation which in turn gives the 
equations of state. 

It is tempting to speculate that this dynamics is related to the renormaliza- 
tion group evolution in quantum field theory: the thermodynamic phase space 
is 'theory' space, the coupling constants being the co-ordinates on it. Partial 
support for this conjecture is provided by the fact that the 'time' variable of 
our dynamics is related to ( but is not always the same as) entropy. 



5.5 The Characteristic Curves of the van der Waals gases 

As an example, consider the case of the van der Waals family: 



F{s,u,v,Pu,Pv) 



vpy - upu + 



P 



(89) 



The characteristic ODEs look quite formidable in these variables. There must 
be a change of variables in which the equations are simple to solve. Put first 



g'' = s + - logw, q^ — V, q'' = uv 
to simplify the formulas. The conjugates are given by comparing 



(90) 



ds-pudu-pydv ^ dq° -pidq^ -P2dq'^ ^ pi = Py + -^--Pu, P2 = — (91) 

ZV V V 

Then (after factoring out an overall v^e^^^^) the hypersurface can be de- 
scribed by the vanishing of 



FiqO,~q\f^p,^p,)^e^?Oplpl-^. 



(92) 



Now we make one more change of variables 



Pi = e^^ Pi, P2 — P2 



(93) 



and choosing q", q , q such that 

a ^ dq° -pidq^ -picJ<?2,=> 9° = 



9^ = <7\ q'^^f, 



(94) 



*The equations in Ref. [3] don't have the term proportional to F. On the surface _F = 
our equations agree with theirs. 
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27 

F=pIpI-y- (95) 

Because F is independent oi q^,q^,q^ in these variables, pi,P2 are conserved 
quantities, and q'^,q^,q'^ depend on 'time' t linearly: these are the 'normal co- 
ordinates' of the characteristic equations. Moreover from the homogeneity of 
-r -r g , 

^=2p.rf. f =3pM (97) 

Transforming back to the original notation, we see that v = q^ and uv = q^ 
evolve linearly in 'time' t.: 

27 

v = q^ + 2-pxvlt. uv = ql^Zv\vlt, 9° = ^y* (98) 
If we identify the constant of integration 

9o = ^ <& = -« (99) 

and eliminate t we get 

uv -\- a 3 pi 3 pi 3 v-py — upu + | 



2p2 2p2 2 pu 



(100) 



Transforming to the original intensive variables,we get one form of the equation 
of state: 

uv + a _ I^-Wt + I 



V — b 



T 

1[vP-u+It] (101) 



Others forms follow similarly. 

What is the meaning of the 'time' variable in this evolution, in terms of the 
original variable? We get (apart from a constant fixing the origin of t) 

t=^^q° = ^vh^'. (102) 
5 X 27^ 27 ^ ' 

The characteristic curves do not change if we replace t by any monotonic func- 
tion of itself. Thus we can think of s + | log?; cqiiivalently as the time variable. 

The thermodynamic 'time' variable is thus related to entropy, but is not 
exactly the same as it. Nevertheless, it seems reasonable to think of it as a 
measure of the coarse-graining of the underlying microscopic system. 
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6 Thermodynamics of a Black-hole 



It is well-known [16] that a black hole has an entropy proportional to the area 
of its horizon. Einstein's equation imply that the total area of horizons cannot 
decrease as long as the source of gravity satisfies the positive energy condition 
(a weak form is sufficient). 

The first indication of an analogy between thermodynamics and black holes 
was found in the work of Christodoulu and Rufhni (17j who were analyzing the 
Penrose process of extracting energy from black holes. Beckenstein [20] 
and Hawking^^ showed that this was more than a mere analogy: black holes 
must have a temperature and an entropy when quantum effects are taken into 
account. For a recent review including a historical prespective, see Ref. [18] . 

Strominger and Vafa [22" have shown how to account for black hole entropy 
from string theory. Other approaches to involving loop quantum gravity [23], 
non-commutative geometry |24| are being vigorously pursued as well. 

For simplicity let us restrict ourselves to spherically symmetric black holes. 
The metric can be brought to the form 

ds^ = VMdt" - - r^dn"". (103) 

V{r) 

The horizon occurs when V{r) = 0. The temperature of this horizon is propor- 
tional to the acceleration of the null Killing vector, which works out to 

K = ir(ro) (104) 

The area of the horizon is, of course, 

A^A-Krl. (105) 

Einstein's equation of motion for the metric imply that any infinitesimal varia- 
tion of these quantities are related by an analogue of the thermodynamic law: 

dM = —ndA (106) 

where M is the mass of the black hole multiplied by Newton's constant. ( In 
units where c = 1, this has dimensions of length.) Thus A is corresponds to 
entropy, temperature to ^ and M to internal energy. This model of a black 
hole has one thermodynamic degree of freedom. 

6.1 The Schwarzschild black hole 

The simplest case, the Schwarzschild black hole, has 

2 M 

V{r) = 1 - — (107) 
r 
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It has the two equations of state 



K=^, A=1&-kM'^ (108) 
4M ^ ' 



Note that 

j^_dM__ (dA_Y^ 



(109) 



as needed. 

These equations of state are analogous to the ideal gas laws we discussed 
earlier. They describe a curve in the three dimensional thermodynamic phase 
space (M. A, k). In order to find a Hamilton- Jacobi theory, we must depart from 
this ideal case by introducing some parameter the deforms Einstein's equation. 



6.2 SchwEirzschild-Anti de Sitter Family 

The choice of such a deformation is a question of physics; the most obvious one 
physically is a cosmological constant. 

If the cosmological constant is positive, even in the absence of mass there 
is a horizon; a black hole will have two horizons, with associated temperatures 
and entropies. Since our aim now is to keep things as simple as possible, we will 
choose the cosmological constant to be negative. 

A = -|. (110) 
In the absence of mass the solution to Einstein's equations 

- ^9iJ.u = i^Qiiu (111) 



is the AdS ( Anti de Sitter) metric 

/V"" 1 + 



ds^ = I 1 + dt" - - r^dn^; (112) 



I.e., 



V{r) = l + -. (113) 

This metric has no horizon as V{r) never vanishes. 

The spherically symmetric solution that has a horizon is 

ds^ = V{r)df r-'dn^ (114) 

V{r) 

where 

2M 

V{r) = I - — + (115) 
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The horizon radius ro is in 1-1 correspondence with the mass: 



2M = ro + 

since they are both positive. The area of the horizon is 



and the temperature is where 



(116) 



(117) 



(118) 



Wc can now eliminate the cosmological constant to get the surface in the ther- 
modynamic phase space describing the Schwarzschild-AdS black holes: 



3M 

47r 



and the Hamilton-Jacobi equation 



3M- 2y4 



dM 
'dA 



(119) 



(120) 



In this family, the Schwarzschild black hole has the largest entropy for a given 
mass. In this sense it is a kind of 'ground-state'. The solution of this differen- 
tial equation depends on one constant of integration, which is the cosmological 
constant. 

6.3 The Chciracteristic Curves of a Black hole 

The above equation give the hamiltonian 



F{M,A,pA) = {iM-2ApA)- 
Again, a Legendre transformation will simplify the dynamics. Put 



A = 47re2«, q° = e'^'^M. 



(121) 



(122) 



Then 



dq° -pdq = Q^p = e'^" [lApA - 3M] . (123) 

After factoring out an overall factor of — e^' the Hamiltonian of the S-AdS black 
hole becomes: 

F{q\q,p)=p + e-^''. (124) 
The dynamical equations are, in these variables, 



^ = 1 * = 2e-29 ^ 
dt ' dt ' dt 



(125) 
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The solution with F = is 



2t 



,0 



1 



2t 



(126) 



q 



— e 
2 



+ a 



where a is a constant of integration. Converting to the original variable M = 
e^'g^ and noting that e'' = ro we get the equation of state 



Thus we identify the constant of integration with the cosniological constant. 



In this picture the 'time' variable of thermodynamic evolution is the logarithm 
of the entropy of the black hole. 

In general a black hole can carry also angular momentum and electric charge. 
The family of such black holes will be described by deformations of Einstein- 
Maxwell equations at short distances. The obvious deformation of Maxwell's 
theory would be the Born-Infeld action. For the Einstein Lagrangian itself 
there are several higher derivative terms that can be added. But the analysis 
gets quite complicated as we have to solve for the black hole metric. It should 
be interesting however. 

Finally, recall that the Hamilton-Jacobi formulation of mechanics gives the 
shortest route to quantum mechanics. We have considered elsewhere the possi- 
bility of a quantum thermodynamics j25j . 
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